Casimir momentum of magneto-chiral matter 
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We consider a scattering formulation of radiative momentum transfer to a Magneto-Chiral system 
that is subjected to a constant background magnetic field. The system takes the form of a collection 
magnetic dipoles that exhibit the Faraday effect. It is shown that the first non-trivial contribution 
of the momentum transfer to the object from the radiation field occurs at fourth order in the Born 
series. 
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Casimir and dispersion energies [l|, [|| as originally de- 
rived in the context of two parallel conducting plates or 
dielectric bodies, are by now a well understood facet of 
quantum electrodynamics with measurable experimental 
implications. The ramifications stretch a broad range of 
science, from the presence of a cosmological constant [|[ , 
to stiction forces in nano-scale architecture and applica- 
tions [4J. See [5[ and [6| for recent surveys. 

In the above scenarios, one usually calculates an energy 
or a force starting from the energy-momentum tensor of 
electromagnetism. Of course the measurable observable, 
that is the energy obtained by integration over all space, 
is but one component of a full tensor of observables. In 
particular, the time-space components are the associated 
momentum density of the fields which upon integration 
over all space gives the total momentum contained in the 
field Q. In the case of a single parallel plate, it is pre- 
cisely because the incident vacuum photons from the left 
and right tend to cancel that the plate does not develop 
a net momentum. This preserves the translational in- 
variance of the vacuum. A natural question then to ask 
is how to find a system whereby a momentum can be 
displayed. 

Magneto-Chiral (MC) systems are bodies that display 
electromagnetic response functions that depend linearly 
on an external magnetic field (54l0j. In addition to this, 
the object is also chiral (breaks parity symmetry). This 
could be because the object as a whole has no reflection 
symmetry, or because the electromagnetic response func- 
tion has spatial dispersion. With much of this structure 
in mind, there have been recent studies [lll - [l7j to see 
if momentum transfer from the vacuum to bodies with 
complex media attributes and subjected to external elec- 
tromagnetic fields is possible. The general conclusion 
seems to be that a momentum can be transferred to a 
body, though the relative magnitudes of the momenta 
and assumptions appear not to be consistent with one 
another. 

The main result of this letter is that we demonstrate for 
a collection of magnetic dipole scatterers that receive a 



classical perturbative correction to its permeability from 
an external magnetic field, a non-zero momentum trans- 
fer to the body as a whole. It is a fourth order pertur- 
bative result and requires at least four scattering centres 
to be present and held in a rigid configuration. Further, 
they should be arranged so that the resulting tetrahe- 
dron (with the four particles placed at the vertices) has 
no parity symmetry so that the vacuum photons get to 
see a chiral structure. If these conditions are fulfilled 
then a non-zero momentum develops which scales as the 
fourteenth inverse power of the length scale of the tetra- 
hedron. Whilst the numerical value of this momentum 
for the pure quantum vacuum case is far too small to be 
measured experimentally, there are hopes of being able 
to measure its classical counterpart. The calculations are 
performed using the notation and conventions in 18[ and 
multiple scattering theory. 
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FIG. 1: The Magneto-Chiral object consists of four point par- 
ticles in a background complex dielectric, that have a mag- 
netic polarisability a M (w) located at positions X x , ■ ■ ■ ,X A . 
The coordinates of the point particles form the vertices of 
a tetrahedron T. This is subjected to a constant external 
magnetic field B°. 



a. Radiative momentum transfer formulation. The 
system that we investigate is illustrated in Figure [1] It 
consists of four magnetic dipoles held rigidly in position 
with respect to one another that is placed in a constant 
external magnetic field (that has been switched on adia- 
batically) and a background dielectric that is a complex 
function of frequency. The question we now try to answer 
is: Will a radiative contribution to the momentum of the 
body result, in particular for the quantum vacuum ? 
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In order to evaluate the quantum electromagnetic field 
momentum for the proposed setup, we need to employ 
a theory that deals with the situation of light-matter in- 
teractions in a consistent way. The framework in which 
we work is that of macroscopic quantum electrodynam- 
ics (see for instance [IH), which ensures in particular 
that the equal-time commutation relations are main- 
tained. One could also envisage performing the calcu- 
lation within a relativistic QFT setting if the potentials 
can be generalised to appropriate scalar fields. The field 
momentum derived from the classical Lorentz force law 
is (at some instant) 



with 



Mv + P(t) = constant, 



P(t) = J d 3 x(E(t, x) A B(t, x)). 



(1) 



(2) 



We now pass to the quantum version of Equation ([2]l 

<P(t)) = J d 3 x(0\E(t,x) AB(t,x)|0), (3) 
with the spectral decomposition 

/>oo 

E(t,x) = dwE(ui,x)e~ iut + H.c, (4) 
Jo 

and similarly for B(t, x). They are linked by the Faraday 
equation 



B(x,u) = V A E(x,u), 



(5) 



together with its hermitian conjugate. By making the ap- 
propriate substitutions and point splitting the arguments 
this can be recast as 



poo />oo /> 

(P a ) = dw dJ \ d 3 x lim (-i/w)((0|Ej(a/,w / )^E 6 (x,w)|0> - (0|Ej(x / ,w')^E a (a;,w)|0» 

JO Jo J x'-tx 



+ c.c. 



(6) 



An important point to appreciate here is that the per- 
meability has both a reciprocal part and a non-reciprocal 
part due to the external magnetic field. Although the 
non-reciprocal parts arise from the point interaction 
where the external magnetic field couples to the scat- 
tering object, this is enough to require the correlation 
functions and thereby the Green tensors, to take a differ- 
ent form. Note here that we are choosing to work with the 
momentum density associated with the canonical energy 
momentum tensor rather than the Poynting vector; the 
latter is expected to integrate to zero [9J so a permeabil- 
ity is required to distinguish the two vectors. Using the 
expectation values of the noise currents at zero temper- 
ature, we find for the two point function of the physical 
fields:- 

(0|E B Or',a/)Ej(*,w)|0) = -^L-8{lo' - lo) 

x(G ab (x',x\u)-G ba (x,x'\uj)% (7) 

with (0|E£(a;', uj')E b (x, w)|0) = 0. Of course for a stan- 
dard reciprocal permittivity, the above simplifies to the 
imaginary part of the Green tensor. 

b. The Green tensor multiple scattering expansion. 
It remains now to evaluate the Green tensor. This can 
be done using a multiple scattering expansion. The full 
Green tensor satisfies the vector Hclmholtz equation 

(V A fi-\x) • V A G{x,y\u)) ab - — (e • G) ab (x, y\w) 

c 2 - 

= 5 ab 5 3 {x-y),(8) 



where we choose e = e ab (x,uj) = S ab for the back- 
ground relative permittivity, whilst the relative perme- 
ability of the MC object n ab {x,bS) must now account 
for the magnetic scatterers. Each scatterer has the form 
[i ab {x,u,)-6 ab = S^x-x 1 )^ 1 {u)8 ab + iuP M {u)e abc B% 
where a M (to) and f3 M (uj) are the magnetic polarisability 
of the atomic system and a correction due to the ex- 
ternal magnetic field, whilst x 1 is the position of the 
J-th particle. The magnetic polarisability is given by 
the expression a M (uj) = a M (0)uq / (w^ - uj 2 + iYuj). By 
use of the Lorentz force law as applied to a bound sys- 
tem of charges that is overall charge neutral, one finds 
the relation f3 M (uj) = a M (u)) 2 /e. A further remark 
about the dimensions of the various quantities is that 
the combination |B°|a M (o;)cL>/e is dimensionless. The 
relative permeability can be turned into the correspond- 
ing T-matrix [8j, together with the multiple scatter- 
ing expansion of the Green tensor (in symbolic form) 
G = G° + G°TG° + G°TG°TG° + • • • , where G° is 
the free vacuum Green tensor. We shall be specifi- 
cally interested in the fourth order term corresponding 
to the scattering from all four particles. An example 
set of scattering events from y to x contributing to the 
fourth order Green tensor is G°(x, X 1 )TG°(X 1 , X 2 )T 
G°(X 2 , X 3 )TG°{X 3 ,X 4 )TG°{X 4 , y). 

As a final step in evaluating Equation (j6|) we Fourier 
transform to momentum space 

{Pa) = (^)i duju IwT 3lka 

x [G bb (k,k\LJ,B°)-G bb (k,k\w,B )*]. (9) 
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A notable feature of the above expression is the absence 
of a transverse contribution. The Green tensor can also 
be split into two distinct pieces. One is the standard 
reciprocal part due to multiple scatterings from the re- 
ciprocal potentials a M and a second piece that is non- 
reciprocal and proportional to B°. The reciprocal piece 
does not contribute to the momentum. 

One is typically interested in evaluating a rotation- 
ally averaged version of this momentum, P a , particu- 
larly from an experimental stand point. A convenient 
way to do so mathematically is to form the scalar prod- 
uct (P a )B° and then extract from this the bilinear term 
B°B° using Equation ©. From here the replacement 
B°B° -> l/3(B°) 2 <5 a6 is then made. The rotationally 
averaged momentum then takes the form 



Pa = gB° a , 



(10) 



where g is a pseudo-scalar. In the following section, we 
will evaluate g = |P Q |/|B°|. One should also compare 
this with the result found in [l2j for chiral molecules in 
a magnetic field. 




H(Hz) 



FIG. 2: A plot of the spectral momentum density (in units 
of kgmT -1 ) as function of imaginary frequency (Hz) for a 
collection of four Sodium atoms with parameters L — lOOnm, 
ojo = 3.1 x 10 15 s-\ 



c. Results and magnitude of the interaction. At this 
point it is necessary to evaluate Equation © either sym- 
bolically or numerically. We have used Mathematica for 
this purpose. Before doing this however we can make 
some analytical statements. Firstly there are no one, two, 
or three body amplitude contributions to Equation ©. 
Essentially this is because after the rotational averaging 
we are left with one three dimensional Levi-Civita tensor 
that has the effect of forming a determinant of the par- 
ticle separation vectors (typically a pseudo-scalar of the 
form X_ 12 ■ (X 23 A A 34 )). It is only at fourth order that 
this can be non-zero i.e. it is a necessary condition that 
three of the separation vectors form a linearly indepen- 
dent set and the tetrahedron T has a finite volume. This 
is of course why we have chosen to consider a system of 
four particles from the outset. A second point is that 



the integral over all space in Equation ([6]) can be done 
resulting in a simple product of free space propagators 
over the edges of T that form a loop. These in particular 
occur as the overall phase factors in the final expression 
as the length of the closed loops due to the compounding 
of the individual propagators. It necessarily results in a 
closed loop of propagators, that start and finish on the 
same scattering site. 

This results in the final form (with k :— uj/c) for the 
pseudo-scalar density 



dg_ 



Im 



4h 

3e^ 



n=l 



2 e ikL, 



F n (ikL) i 



,M 



' 11 



{kc)n f 



(11) 



where L n is the length of one of the three inequivalent 
closed loops (e.g. scatterings from particles in the order 
12341 with n indexing the loop). In addition, L is the 
characteristic length scale of T (taken to be the length of 
the shortest edge), and F n (ikL) is a dimensionless poly- 
nomial of degree < 8. We are neglecting recurrent scat- 
tering and working to fourth order. Then the pseudo- 
scalar is given by g — J °° duj(dg / du). Equation ([TT|) can 
be explicitly evaluated for a collection of different shaped 
tetrahedra. Notably it vanishes when T is chosen to be 
regular i.e. parity invariant. At low frequency 



dg_ 

did 



AH 



3efioL 



li 



(c^((%o) 5 (^) 2 . (12) 



An example evaluation for a set of vectors X_ l = 
(0,0,0), X 2 = (£,0,0), X 3 = (0,2i,0), 2U = 
(L/2, L/2,3L) gives a non-zero result. In Figure [2] we 
give a sample plot for Sodium atoms of the spectral mo- 
mentum density as a function of frequency having per- 
formed a Wick rotation of frequency to the imaginary 
axis. 

The frequency integral can be evaluated in the retarded 
limit (where we replace the polarisabilitics with their 
static values) resulting in 



fJ 



( Ahc\ 
\3enoJ 



, M 



(0)W) 



1.4 x 10 4 
7Ji 



(13) 



which with the parameters from Figure [5] and |B° = 
10T give a momentum \P a \ « 3.1 x 10~ 47 kgms^ 1 with a 
corresponding MC body speed of v « 2.1 x 10~ 22 ms -1 . 

It is also possible to evaluate the momentum at fi- 
nite temperature using the standard Matsubara series 
by inserting a factor of coth(hu / (2k bT)) into Equa- 
tion (ITT|) . The Matsubara frequencies are given by fl m — 
2TtkBTm/h where m is a positive integer. In the low tem- 
perature limit this can be evaluated for the first few Mat- 
subara frequencies. This is done in Figure [3] for temper- 
atures less than 300K. An order of magnitude estimate 
is also straightforward: the ratio of (g(T) — g(0))/g(0) is 
approximately for the first Matsubara frequency 



g(r)- g (o) 

9(0) 



a M ^i)\ 5 (Lk B T^ 3 



a M (0) 



\ he 



(14) 
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g(T)-g(T=0) 
(xlO- 50 kgm/(sT)) 




FIG. 3: A plot of the momentum difference per Tesla as a 
function of temperature for a collection of four Sodium atoms 
with parameters as in Figure 2. This is the low tempera- 
ture regime determined from the length scale L for the first 
Matsubara frequency. 



where N\ is a numerical constant arising due to the 
different frequency integration being performed at zero 
and finite temperature. Using the data for sodium gives 
fa 0.5% correction at T = 300K. 

The high temperature limit can also be calculated 
which is defined by T ^ hcKksL). In this case the ther- 
mal factor at imaginary frequency is cot(hQ/2kBT) w 
IksT jfr£l. The frequency integral can be done as before 
with one less power of frequency resulting in 

lim 9(T) « (^%) (a M (0)Mo) 5 , (15) 

where the S's have cancelled out and we enter a regime 
of classical linear temperature dependence. 

So far the calculation we have performed has been for 



a quantum system of noise currents that have vacuum 
expectation values proportional to h. One can also form 
a classical counterpart of Equation (fTTj) by considering a 
classical source. For an isotropic distribution of classical 
radiation such as from a laser, (dgfduj)dui must be re- 
placed by Ag = dg / 'du)(F/ c)(2fiw 3 / (nc 3 ))" 1 , where F is 
the power per unit area of the laser source; the latter fac- 
tor is the energy density per unit frequency of the vacuum 
photons we are replacing. An order of magnitude esti- 
mate for the momentum can be obtained by evaluating 
the expression off- resonance (at uj = l.OOlwo) to avoid re- 
current scattering effects not incorporated here. We use 
\P a \ « |B°|(A ff | w=1 .ooi Wo . With values of F= l&Wm~ 2 
and B°| = 10T, this is found to be \P a \ w IQ-^kgrns' 1 
with a corresponding speed v w 10 ms . The effect 
can therefore be greatly enhanced by the use of classical 
sources as compared to the vacuum case. 

d. Summary and perspective. In this letter we have 
investigated if a simple Magneto-Chiral object can de- 
velop a momentum from the quantum vacuum. This 
effect is allowed by symmetry. In itself, the result is a 
tiny effect for the quantum case, but nevertheless non- 
zero for the case of magnetic dipoles. For electric dipolcs 
it vanishes rigourously. Thermal corrections have been 
calculated in the low temperature limit, as well as a high 
temperature limit. The real frequency form has also been 
investigated in a classical setting and the relative en- 
hancement of momentum found off-resonance due to a 
classical driving source. 
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